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Analyzing recent ALICE data on inelastic pp scattering at the LHC energies we show that charged 
particle distributions exhibit geometrical scaling (GS). We show also that the inelastic cross-section 
is scaling as well and that in this case the quality of GS is better than for multiplicities. Moreover, 
exponent A characterizing the saturation scale is for the cross-section scaling compatible with the one 
found in deep inelastic ep scattering at HERA. Next, by parametrizing charged particles distributions 
by the Tsallis-like formula, we find a somewhat unexpected solution that still exhibits GS, but differs 
from the ’’standard” one where the Tsallis temperature is proportional to the saturation scale. 


I. INTRODUCTION 

It is believed that gluon distribution inside a hadron 
saturates at small Bjorken x (see Refs. mm for review). 
This is a consequence of the non-linear QCD evolution 
equations, known as Balitsky-Kovchegov (BK) equation 
[3j or in a more general case as JIMWLK equation [1], 
that poses traveling wave solutions [5]. The latter prop¬ 
erty, in the QCD context, is called geometrical scaling 
(GS) (6]. An effective theory relevant for the small 
Bjorken x region is so called Color Glass Condensate 
(CGC) [7]. For the purpose of present work the details 
of the saturation are not of primary importance; it is 
the very existence of the saturation scale which plays the 
crucial role 

Geometrical scaling means that some observable N 
that in principle depends on two independent kinemat- 
ical variables, say x and Q 2 , in fact depends only on a 
specific combination of them denoted as r: 


N(x,Q 2 )= JF(t). 

(1) 

Here function T in Eq. (JT|) is a dimensionless function of 
a scaling variable |8] 

r = Q 2 /Q 2 (x). 

(2) 

and 


Qs( x ) = Qo ( x / x o)~ x 

( 3 ) 


is the saturation scale. The power law form of the satu¬ 
ration scale is dictated by a saddle-point solution to the 
BK equation and has been tested phenomenologically for 
different high energy processes B-m 

Here we are coming back to pp scattering [9] in the con¬ 
text of recently published ALICE data m for charged 
particle distributions at three LHC energies 0.9, 2.76 
and 7 TeV. After discussing shortly in Sect. [II] how GS 
emerges in the h? factorization scheme, we shall show in 
Sect. |III| that recent ALICE data indeed exhibit geometri¬ 
cal scaling with, however, A exponent being different than 
in the case of deep inelastic (DIS) ep scattering. Interest¬ 
ingly, we shall also show that the inclusive cross-sections 
scale somewhat better and with an exponent that is very 


close to the DIS value A = 0.32 [TO]. This result calls for 
better understanding of the impact parameter picture of 
pp scattering in the context of the saturation physics and 
the Color Glass Condensate theory. 

Another topic addressed in the present paper is the 
shape of the scaling function introduced schematically in 
Eq. |l]). Function T can be in fact obtained only nu¬ 
merically within some specific model. Here, we shall use 
phenomenological parametrization in the form of Tsallis- 
like distribution M applied successfully in the past to 


we briefly describe how GS should be reflected in the 
Tsallis distribution. Next, in Sect. [V] we shall try do 
fit Tsallis-like parametrization to the ALICE data. Un¬ 
fortunately, as already remarked in the original ALICE 
publication m , this piece of data does not admit good 
quality Tsallis fit. Nevertheless, we invoke a procedure 
that allows for rather good description of the data in 
the range of moderate transverse momenta where GS is 
expected to occur. Somewhat unexpectedly we find GS 
scaling solution that is very different from the ’’standard” 
one described in Sect.|IV| Unfortunately GS in this solu¬ 
tion is rather accidental and will disappear at very high 
energies. Whether this is only a property of the Tsal¬ 
lis parametrization ’’forced” to describe ALICE data, or 
a real prediction, remains to bee seen. We conclude in 
Sect. EJ 


describe spectra of charged particles [14] —[16] . In Sect IV 


II. GEOMETRICAL SCALING IN HADRONIC 
COLLISIONS 


The cross-section for producing a moderate px gluon 
in hadronic collision can be described as El¬ 


der 


C 


= — drR t a s (fc4)<pi(^i, k^)ip 2 {x 2 , (k - p )x) 
dyd z p T PtJ 

( 4 ) 


where C contains color factors and numerical constants. 
Bjorken x’s of colliding partons read 


x l,2 


Pr ±y 


( 5 ) 
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however, since in the following we will be interested 
in central rapidity production, i.e. y — 0, we have 
X\ ~ X 2 = x. Functions ipi t 2 are the unintegrated gluon 
distributions in hadron h\ and /12 respectively. There 
exist many models of unintegrated gluon distributions 
(see e.g. d8j); the most simple ones are the one of the 
Golec-Biernat-Wiisthoff (GBW) model [TS] or the one 
by Kharzeev, Levin and Nardi }20| . They share two im¬ 
portant features: geometrical scaling and dependence on 
the transverse area S± whose precise meaning is best un¬ 
derstood in a picture where also the impact parameter is 
taken into account EQ OH]. Therfore 


note that GS is in fact a property of Eq. (J7| and that the 
multiplicity scaling (101 is based on |9]) which is not so 
obvious for the scattering of small systems, like pp. 

In order to integrate (10) over d 2 pT we have to change 
variables 


PT = ( 11 ) 


where the average saturation scale is defined as 


Q S (W) 


Q 0 


/£IF\ A / (2+A ) 


( 12 ) 


<p(x, £f) = S^(j){k\/Q 2 s (x)) (6) 


where <fi is a dimensionless function of the scaling vari¬ 
able k^/Q 2 {x), rather than independently of x and k^,. 
Of course geometrical scaling is only an approximation 
and is expected to break for large Bjorken a:’s and also 
for large transverse momenta. We also expect GS break¬ 
ing for small fcx where non-perturbative effects including 
effects from the pion mass are of importance. Ignoring 
these effects and neglecting also momentum dependence 
of the strong coupling constant we arrive at 


da 

dyd 2 p T 



( 7 ) 


where J- is a universal, energy independent function of 
the scaling variable r: 


Pt = (Pt \ 2 ( Pt \ 

Qsi x ) \QoJ U wj 


( 8 ) 


Here we have used (|3j) for the saturation scale Q 2 . We 
take for Xq = £ x 10 -3 . This implies that in (I8| px is 
in GeV/c and W in TeV. Furthermore for Qq we can 
take without any loss of generality Qq = 1 GeV/c. One 
typically assumes that S± is an energy independent con¬ 
stant. This is true in the case of the GBW model m 
where S± = ao with a 0 characterizing the asymptotics 
of the dipole-proton cross-section for large dipole sizes. 
In the case of heavy ion collisions for fixed centrality, S± 
has geometrical interpretation as an overlap area of two 
colliding nuclei [20]. In this case one can also assume 
that 


da ^ d 2 N 
dyd 2 pT dyd 2 pT 


( 9 ) 


where TV is a multiplicity of produced gluons. Neglect¬ 
ing possible energy dependence of gluon fragmentation 
into hadrons [ 25 ], i.e. adopting parton-hadron duality 
hypothesis [24] , we arrive at: 


1 d 2 N c h 
Pt dydpr 


S±T(t). 


( 10 ) 


Expression ( fl0| ) will be used in the following to look for 
GS in the multiplicity distributions. Let us, however, 


Note that the effective power describing the rise of the 
average saturation scale with energy A e s = A/(2 + A) is 
slightly smaller than A/2. Then 


PrdpT = x-^t-Qs (W)t 2/(2+a) —,. (13) 

2 T A T 


Hence 

dN ch 

dy 


2 + A„ 


Hr) 


r 2/(2+A) 


dr 


x S ± Q 2 s (W) 


= b X SxQ s 2 (W 0 ■ 


(14) 


Here b is an energy independent constant related to the 
the integral of Equation (14) is often used as a 

definition of the saturation scale (with TV c h replaced by 
Ngi U ons) understood as gluon number density per trans¬ 
verse area. 


III. GEOMETRICAL SCALING OF ALICE 
DATA 

In this Section we are going to check whether ALICE 
data |12| on inelastic multiplicity distributions of charged 
particles exhibits geometrical scaling and for what value 
of A. We shall show that indeed GS is reached for A ~ 
0.22 — 0.24 as it is illustrated in Fig. |T] 

In order to find the best value of A we have adopted 
the method of ratios described in more detail in Refs. m- 
Let us denote for simplicity 


N{pr,W) 


1 d 2 N ch 


2ttpt dydpT 


w 


(15) 


We form ratios of spectra expressed in terms of the scal¬ 
ing variable r rather than in terms of pt' 


R 


w/w 


( T ) 


N(t,W) 
7V(t, W') 


(16) 


and request that R ~ 1 over the largest possible interval 
of r. Note that this method is sensitive only to the value 
of A and not to the actual values of parameters Qq and xq 
In the present case we choose 7 TeV for W and W\ = 2.76 
or W 2 =0.9 TeV for W'. Therefore we can form two such 
ratios, which are depicted in Fig. [ 2 ] for A = 0 {i.e. for 


















3 


4 = 0 /l = 0.24 




FIG. 1. Charged particle spectra measured by ALICE m , plotted as functions of pt (left panel) and as functions of the scaling 
variable r Q for A = 0.24 (right panel). Black full dots correspond to IF = 7 TeV, red down-triangles to 2.76 TeV and blue 
up-triangles to 0.9 TeV. Solid lines correspond to the Tsallis fits from Sect. [V] 
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FIG. 2. Ratios of charged particle spectra measured by ALICE ca, plotted as functions of pt (left panel) and as functions 
of the scaling variable r © for A = 0.24 (right panel). Red down-triangles correspond to to the ratio 7/2.76 TeV and blue 
up-triangles to 7/0.9 TeV. 






y/r = pt) and for A = 0.24. We see that indeed, these 
two ratios that rise rather steeply with pt , remain flat 
and close to 1 if plotted in terms of y/r for A = 0.24. We 
interpret this as a signature of geometrical scaling. 

In order to decide on the best value of exponent A we 
need to provide a quantitative criterion measuring the 
’’average distance” of experimental values of Rw/w from 
unity. To this end we propose the following procedure. 
Since for A values relevant for the present analysis the 
first few R points corresponding to low pt lie above 1 
(which is the sign of GS violation in a region when non- 
perturbative effects are of importance) we pick up the 


first point for which 

Rw/Wi ,2 ("Tstart) If/ A // ( 7"start) ■ 

Here A r is an experimental error of R. For points with 
t > T s tart ratio R is either close to 1 within the experi¬ 
mental errors, or it is falling below 1 exceeding A#. Next, 
since for large transverse momenta pt spectra are getting 
harder with increasing energy, the values of R start to in¬ 
crease with t, getting again larger than 1. This is well 
visible in Fig. [3] where the vertical scale has been magni¬ 
fied with respect to Fig. [2] for better resolution. 

Starting from r sta rt that of course depends on energy 
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FIG. 3. Same as Fig. [2] with different scale for better resolution. Solid lines in the right panel correspond to the Tsallis fits 
from Sect. [V] 


W-[ . 2 , we compute mean square deviation for given A 


^Wi, 2 W 


nwi. 


E 


, (A) ^ 

“ V ' T n =T s atrt 


(R(r n ,W h2 ) ^ if 
A|j(t„, Wl j2 ) 


(17) 


where ^^(A) is a number of points between T sta rt 
and r en( j. We increase r en d up to the last point where 
^Wi 2 W < 1- In this way we obtain n^(A) which is the 
number of points that exhibit GS for given W' = W\ ^ 
and for given A, which are plotted in Fig. [4] Now we 
look for maximum of nwi(A) + riw 2 (A). This happens 
for A = 0.24. As seen from Fig. [3] W 2 = 0.9 TeV points 
scale in a shorter interval of r, which translated back to 
transverse momenta corresponds to p™ ax = 3.1 GeV/c. 
We see from Fig. [ 3 ] (right panel) that although Syy 2 < 1 
there are 0.9 TeV points (blue up-triangles) in the region 
Tstart < t < T en d which are below 1 outside the experi¬ 
mental error. If we demand that all points between T sta rt 
and r en d should be equal to unity within experimental 
errors, then A = 0.22. This is the value of A used in 
Refs. [3S] and the relevant plot is shown in the left panel 
of Fig. [3] The corresponding p™ ax is shifted down to 
2.9 GeV/c. 

Interestingly, when we repeat this procedure for the 
cross-sections which are obtained by multiplying the 
multiplicity spectra by the minimum bias cross-section 
ctmb(W) given explicitly in Ref. [12], we find that GS oc¬ 
curs at a higher value of A = 0.31 — 0.33. This by itself is 
not surprising since <7 mb(W) depends on energy and this 
dependence makes A different than in the case of multi¬ 
plicity. What is, however, surprising and encouraging is 
that the value of A is now consistent with DIS. Moreover, 
the range of GS is now larger, up to p™ ax = 4.25 GeV/c. 
This is depicted in Fig. [5] The explanation of this obser¬ 
vation is beyond the scope of the present paper, however, 
it is clear that it requires a more sophisticated model of 
S± of Eq. ([ 9 ]) , which in the present analysis is assumed to 
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FIG. 4. Number of points that contribute to (171 for two dif¬ 


ferent ratios: 7/2.76 TeV (red down-triangles) and 7/0.9 TeV 
(blue up-triangles) plotted as functions of A. 


be an energy independent constant in the case of multi¬ 
plicity scaling or minimum bias cross-section in the case 
of cross-section scaling. 

We have performed similar analysis for the UAl 
data |2S| for pp cross-section at yfs = 0.2, 0.5 and 0.9 TeV 
with similar result that A « 0.34. Here, however, the 
data extends only up to ~ 7 GeV/c (for two lower ener¬ 
gies) and the tail is quite noisy, namely the ratios of the 
cross-sections fluctuate quite significantly for yfr > 3. 


IV. GEOMETRICAL SCALING AND TSALLIS 
PARAMETRIZATION 

It is well known that particle spectra at low and 
medium transverse momenta can be described by ther¬ 
mal distributions in transverse mass m t = \/p\ + m 2 
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FIG. 5. Ratios of charged particle cross-sections measured by ALICE, plotted as functions of the scaling variable r |8| for 
A = 0.31 (left panel) and for A = 0.33 (right panel). Red down-triangles correspond to to the ratio 7/2.76 TeV and blue 
up-triangles to 7/0.9 TeV. 


with ’’temperature” T which is a function of the scat¬ 
tering energy [27] . It is also known that more accurate 
fits are obtained by means of Tsallis-like parametrization 
m where particle multiplicity distribution takes the fol¬ 
lowing form (see e.g. I2H1 ): 


1 d 2 N ch dN ch p C n 


27rpx dydpT dy E 2ir 


1 + 


Et 

nT 


(18) 


where Et = \J to 2 + Pt — m. In what follows we shall 
keep to = 0 what implies p/E = 1. Here 


C n = 


(n — 1 )(n — 2) 
n 2 T 2 ' 


(19) 


of ([18]). Indeed 
dN ch 


Coefficient C n in Eq. (|19|) ensures proper normalization 

d 2 N ch 


dy 


/ — 

J 2np T dydpT 

OO 

dN ch 


d pt 


dy 


C n / dp T p T 


1 + 


Pt 

nT\ 


( 20 ) 


where the last integral is equal 1 /C n . Here n and T are 
free fit parameters that depend on particle species and 
on energy. 

In the limit n —> oo (or equivalently for small pt) dis¬ 
tribution (18) tends to the exponent 


1 d 2 N ch 
Pt dydpT 


dN ch 1 
dy T 2 


exp(—p T /T). 


( 21 ) 


Substituting (111 into (21) we arraive at: 

1 d 2 N ch ^ dN ch 1 / Q S (W) 

Pt dydpT dy T 2 (W) ^ \ T(W) 


■1/(2+A) 


Equation (22) exhibits geometrical scaling exactly, only 
when [29] 


TiW) = —Q S (W). 

K 


Then, using Eq. (14), we get: 


1 d 2 N, 


ch 


bx 2 


Pt dydpT Q\ 


exp . 


(23) 


(24) 


Indeed, (24) is energy independent. This would gener¬ 


alize to the full Tsallis distribution if exponent n were 
constant. We know, however, from the phenomenologi¬ 
cal fits that n is decreasing with energy making pt spec¬ 
tra harder and - in the same time introducing explicit 
violation of geometrical scaling for particle spectra. 

Let us observe that we can include factor £ into a defi¬ 
nition of b and x, so without any loss of generality we can 
set £ = 1. Therefore we finally arrive at the GS-Tsallis 
parametrization of the pt spectra that reads: 


1 d 2 N ch 


2npT dydpT 


w 


_ B c 
~ n2 n VH 

Vo 




1 + -= 


n\v Q S (W) 


(25) 


where we have introduced new constant B = bx 2 /2 tt and 
explicitly indicated that n is a function of W. This de¬ 
pendence would be of course a source of GS violation. 
Constants B , x and Q 0 should remain energy indepen¬ 
dent. 


V. TSALLIS PARAMETRIZATION OF ALICE 
DATA 

In this Section we are going to check whether one can 
fit ALICE data [T2] with the help of formula (25). In the 


(22) original ALICE paper jT2] it is said that the multiplicity 
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data can be fitted with the Hagedorn distribution m, 
rather than with the Tsallis one. Therefore we could 
expect that ordinary fitting procedures would not give 
a reasonable parametrization of the data. In order to 
enforce Tsallis parametrization we have proceeded in the 
following way. For each LHC energy we have chosen two 
values of px, one in the small px region and one in the 
tail that are displayed in Table hi For pij? w we have chosen 
approximately 0.5 GeV/c that is already above the non- 
perturbative region. For Py Sh we have chosen values that 
are rather far from the end of the spectrum, but already 
large enough to be in the perturbaitve regime. Of course 
our fit parameters do depend on this choice, however, as 
we shall see below, the quality of the Tsallis fits with the 
values of limiting px given in Table [T] is good enough that 
manipulating with these values has not been necessary. 
Let us also remark at this point that our aim here was 
to show certain properties of the Tsallis fits enforced on 
ALICE data at low and moderate transverse momenta, 
since we new from the beginning that this particular piece 
of data does not admit Tsallis parametrization in the 
whole pt range. 


W [TeV] 

Pt w [GeV/c] 

p“ sh [GeV/c] 

0.90 

0.525 

8.5 

2.76 

0.525 

10.5 

7.00 

0.525 

13.5 


TABLE I. Values of px w and p^ Igh used to fit Tsallis 
parametrization to ALICE data (see the beginning of Sect.[V|. 


W [TeV] 

xw 

nw 

B w 

all px 

2 

a W 

Pt < 20 GeV/c 

0.90 

7.0 

8.32 

28.88 

0.46 

0.46 

2.76 

7.9 

7.33 

36.96 

0.87 

0.71 

7.00 

8.6 

6.79 

43.82 

1.49 

0.75 


TABLE II. Parameters entering Tsallis parametrization ( |25| ) 
coming from the fit to ALICE spectra and the corresponding 
values of mean square deviation for all measured px values 
and for px < 20 GeV/c. 



FIG. 6. Mean square deviations defined in Eq.(27l as func¬ 
tions of x. Black circles correspond to W = 7 TeV, red down- 
triangles to 2.76 TeV and blue up-tangles to 0.9 TeV. 


For the px values given in Table [T] we have calculated 
ratios N(p l ™, W)/N(p^ sh , W), both for the data and for 
parametrization (251. In this way normalization param¬ 


eter B canceled out. Now, for fixed value of x we have 
calculated % from the following condition: 


N(p l ™,W) 


N(p^ sb ,W) 


th 


N(p l ° w ,W) 

nW\w) 


n\y. (26) 


exp 


Note that the value of parameter A entering the definition 
of the saturation scale (121 has been already fixed by the 
method desribed in Sect. m Here we have used A = 0.24. 

Next, for each pair (x, %) we have computed mean 
quadratic deviation 


’w 




ir,W)L- N(jfr,w) 


exp 


A 2 (pi r , w) 


(27) 

where i runs over experimental data points at energy W 
up to the maximal px ■ A denotes the experimental error 
of N. The result is plotted in Fig.[6j We see that func¬ 
tions cr(x(x) exhibit minima at three distinct values of 
parameter x. This is the first signal that one cannot fit 
ALICE data with Tsallis distributions that correspond to 
the energy independent parameter x. Therefore we have 


to allow for energy dependent x —> xw which takes us 
away from the geometrically scaling Tsallis parametriza¬ 
tion of Eq.@. 

One can also see that minima of a^y grow with en¬ 
ergy. This is due to the fact, that Tsallis distribu¬ 
tions used here are not able to describe both low px 
part and the very high end of the spectrum simultane¬ 
ously. We have checked that confining the sums in (27) 
to p“ ax ~ 20 GeV/c corresponding to i mSLX = 54 for all 
energies in question, brings down erj^ below 0.8 (see Ta¬ 
ble 0. This means that Tsallis parameterization with 
energy dependent xyy used here is able to describe the 
px spectra at small and moderate transverse momenta, 
be. precisely in the region we are interested in. 

In Table [TT| we collect values of the parameters xw, nw 
and Bw at the minima of a^y. B\y is calculated from 
the condition A7(p^? w , W) | th = N(p l ™, W) | . The re¬ 

sulting spectra together with ALICE data are shown in 
Fig.[l]both for distributions expressed in terms of px and 
in terms of scaling variable \fr. The quality of this fit 
can be also appreciated by looking at the right panel of 
Fig-! where the multiplicity ratios are well reproduced 
without any adjustment of fit parameters. 


Given the fact that the saturation momentum scales 
as a power of energy Q S (W) = (W/Qo) Aeff we have tried 
to fit energy dependence of parameters K\y. nw and % 
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with generic form ao(W/Qo) a with the following result: 

K W = 7.097 (WVQo) 0 ' 1000 , 
n w = 8.199 (W/Qo)-°- w05 , 

B w = 29.76 (W/Q 0 ) 0 - 2013 . (28) 

This result is surprisingly in line with the effective ex¬ 
ponent of the saturation scale which for A = 0.24 reads 
A e ff = 0.1071. Note also that B ~ n 2 (see definition of B 


below Eq. (251), and this dependence is reproduced by 
the fits of Eqs. (281. Although this energy dependence 


follows only from the fit to data, and we do not have any 
model to explain their values, it is a reasonable assump¬ 
tion to take: 


Kw = «o 


nw = no = 


Qs(w) 

Qo 
Qo 


Bw = Bi 


Q s (wy 

Ql(w) 

Qo 


(29) 


where the coefficients ko, no and Bq can be read off from 
Eq. (281. In what follows we shall drop Q 0 = 1 GeV/c 


which was included in (291. Therefore we have: 


d 2 N, 


ch 


dyd 2 p T 


7 * 2 , Tm (n 0 - Q s (W))(n 0 - 2Q S (W)) 


= B 0 Q S (W) 


w 


1 + — Q s 2 (IE)r 1/(2+A) 
n 0 


‘o 

1 -no/Qs(W) 


(30) 


For geometrical scaling to be present we need this 
function to be inedependent of W, i.e. independent of 
Q S (W). For the energies in question (from a few hun¬ 
dreds GeV up to a few TeV) Q S {W) changes from 0.9 to 

1.5. Therefore the factor involving no is in fact close to 1 

_2 

and the main energy dependence comes from Q s (W) in 
front and from the factor in a square bracket in Eq. (301. 


In order to see how GS is reached by Eq. (301 we plot 
in Fig. [7] ratio N(y/r, 7)/N(y/r, W) as a function of y/r 
for W = 0.9, 2.76 and 14 TeV (left panel). Horizontal 
dashed lines at 1 ± 0.15 show 15% band around unity 
which roughly corresponds to the size of the experimental 
errors A R (10 %) and the accuracy of the fit (5 %) - see 
Fig. m GS is present if theoretical solid curves fall within 
this interval. We can conclude from Fig. [7] that with this 
accuracy GS should be seen in the data up to sjr ~ 4 for 
the whole LHC energy range up to 14 TeV. Should Tsallis 
parametrization (30) hold for higher energies we might 


expect shrinking of the maximal yfr where GS is still 
present with increasing energy. Given the fact that for 
fixed r transverse momentum p t is an increasing function 
of energy (see Eq. 0), this may not immediately mean 
that the pt window for GS would be shrinking as well. 

T he s ame conclusion can be reached by looking at 
Eq. (301 where N(y/r, W)/N(y/r, 7) is plotted as a func¬ 
tion of energy for fixed r. In the right panel of Fig. [7] we 


plot N(y/r, W)/N(y/r,7) as a function of W for different 
values of yfr = 0.2, 0.5, 1, 2, 3, 4 and 5 shown next to 
the horizontal axis. 


VI. CONCLUSIONS 

In this paper we have addressed three questions con¬ 
cerning saturation in high energy pp scattering. To this 
end we have used recent ALICE data on inelastic scat¬ 
tering at the LHC m- 

The first question concerned the very existence of ge¬ 
ometrical scaling in multiplicity distributions. By apply¬ 
ing a model-independent method of ratios we have es¬ 
tablished that GS is indeed present in multiplicity spec¬ 
tra over a limited transverse momentum range up to 
~ 3 GeV/c with characteristic exponent A ~ 0.22 — 0.24. 
This exponent is significantly different than in DIS, where 
A = 0.32, and also lower than the one extracted from the 
CMS non-single diffractive data: A = 0.27. We have pro¬ 
posed the solution to this discrepancy by looking at GS 
for the inelastic cross-section rather than for the multi¬ 
plicity distribution. Motivation for this comes from the 
fcx factorized form of the gluon production in pp col¬ 
lisions Q that leads straightforwardly to Eq. (Jt|) and 
from the fact that the proportionality factor between 
the multiplicity and the cross-section © is not energy- 
independent. We have found that the inelastic cross- 
section scales better than multiplicity up to he maximal 
transverse momentum that is larger than 4 GeV/c and 
with the characteristic exponent A ~ 0.31 — 0.33. We 
have also looked at the UAl data for pp scattering and 
obtained similar value of A. We believe that this ob¬ 
servation provides a solution to the discrepancy between 
scaling properties in DIS and in hadronic collisions. 

The second question concerned the universal shape of 
GS and its connection to the Tsallis distribution. We 
have confirmed that the natural answer to this question is 
provided by a parametrization where the Tsallis ’’temper¬ 
ature” T is proportional to the average saturation scale 
Q s (12) and the remaining Tsallis parameter n should 


be an energy independent constant. In practice n does 
depend on energy and this leads to the violation of GS 
for this particular parametrization. 

Finally the third question was whether such a simple 
solution is admitted by the experimental data. We have 
found that recent ALICE data on inelastic charged par¬ 
ticle multiplicity does not admit the above solution, in 
agreement with the original claim of Ref. jT2]. We have 
found another parametrization where Tsallis parameter 
n is inversely proportional to Q s . This parametrization 
indeed exhibits GS in the limited energy range, however 
GS is not obviously extended to higher energies. We have 
concluded at this point that the solution we found was 
rather accidental. It will be therefore interesting to see 
whether this solution will be still present at higher ener¬ 
gies of the LHC run II. 




















FIG. 7. Left panel: Ratios (161 in terms of parametrization of Eq. (301 for W' = 0.9 (blue), 2.76 (red) and 14 TeV (brown). 
Right panel: multiplicities N(y/ri , W) normalized to N at 7 TeV plotted as functions of W in TeV for fixed values of y/ri 
displayed next to the vertical axis. 
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